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The  synapse  number  and  the  related  dendritic  spine  number  in  the  cerebral  cortex  of  primates  shows  a 
rapid  increase  after  birth.  Depending  on  the  brain  region  and  species,  the  number  of  synapses  reaches  a 
peak  before  adulthood,  and  pruning  takes  place  after  this  peak  ( overshoot-type  synaptic  formation). 
Human  mental  disorders,  such  as  autism  and  schizophrenia,  are  hypothesized  to  be  a  result  of  either 
too  weak  or  excessive  pruning  after  the  peak  is  reached.  Thus,  it  is  important  to  study  the  molecular 
mechanisms  underlying  overshoot-type  synaptic  formation,  particularly  the  pruning  phase. 

To  examine  the  molecular  mechanisms,  we  used  common  marmosets  (Callithrix  jacchus).  Microarray 
analysis  of  the  marmoset  cortex  was  performed  in  the  ventrolateral  prefrontal,  inferior  temporal,  and 
primary  visual  cortices,  where  changes  in  the  number  of  dendritic  spines  have  been  observed.  The  spine 
number  of  all  the  brain  regions  above  showed  a  peak  at  3  months  (3  M)  after  birth  and  gradually 
decreased  (e.g.,  at  6  M  and  in  adults).  In  this  study,  we  focused  on  genes  that  showed  differential  expres¬ 
sion  between  ages  of  3  M  and  6  M  and  on  the  differences  whose  fold  change  (FC)  was  greater  than  1.2. 
The  selected  genes  were  subjected  to  canonical  pathway  analysis,  and  in  this  study,  we  describe  axon 
guidance  signaling,  which  had  high  plausibility.  The  results  showed  a  large  number  of  genes  belonging 
to  subsystems  within  the  axon  guidance  signaling  pathway,  macrophages/immune  system,  glutamate 
system,  and  others.  We  divided  the  data  and  discussion  of  these  results  into  2  papers,  and  this  is  the  first 
paper,  which  deals  with  the  axon  guidance  signaling  and  macrophage/immune  system.  Other  systems 
will  be  described  in  the  next  paper.  Many  components  of  subsystems  within  the  axon  guidance  signaling 
underwent  changes  in  gene  expression  from  3  M  to  6  M  so  that  the  synapse/dendritic  spine  number 
would  decrease  at  6  M.  Thus,  axon  guidance  signaling  probably  contributes  to  the  decrease  in  synapse/ 
dendritic  spine  number  at  6  M,  the  phenomenon  that  fits  the  overshoot-type  synaptic  formation  in 
primates. 

Microglial  activity  (evaluated  by  quantifying  AIF1  expression)  and  gene  expression  of  molecules  that 
modulate  microglia,  decreased  at  6  M,  just  like  the  synapse/dendritic  spine  number.  Thus,  although 
microglial  activity  is  believed  to  be  related  to  phagocytosis  of  synapses/dendritic  spines,  microglial 
activity  alone  cannot  explain  how  pruning  was  accelerated  in  the  pruning  phase.  On  the  other  hand, 
expression  of  molecules  that  tag  synapses/dendritic  spines  as  a  target  of  phagocytosis  by  microglia 
(e.g.,  complement  components)  increased  at  6  M,  suggesting  that  these  tagging  proteins  may  be  involved 
in  the  acceleration  of  pruning  during  the  pruning  phase. 
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1.  Introduction 

Primates  show  a  rapid  increase  in  the  synapse  and/or  dendritic 
spine  number  in  infancy  and  early  childhood,  reaching  a  peak.  After 
this  peak  is  reached,  a  decrease  is  caused  by  pruning  [1-5],  This  pri¬ 
mate  type  of  synaptic  change  is  called  overshoot-type  synaptic  for¬ 
mation  [5],  which  is  not  detectable  in  rodents.  Recently, 
abnormalities  in  the  generation  and/or  pruning  of  dendritic  spines 
were  implicated  in  several  mental  disorders  including  autism  and 
schizophrenia.  An  autistic  brain  tends  to  show  more  synapse  forma¬ 
tion  at  the  early  age,  with  less  pruning  after  the  peak  of  the  synaptic 
number  is  reached,  whereas  a  schizophrenic  brain  tends  to  show 
more  pruning  on  an  average,  after  the  peak  of  the  synaptic  number 
is  reached  [6].  To  understand  the  mechanism  underlying  these 
mental  disorders  and  to  find  possible  therapeutic  targets,  it  would 
be  useful  to  identify  molecular  mechanisms  underlying  the  over- 
shoot-type  synaptic  formation,  particularly  of  pruning  after  the  peak 
is  reached.  Pruning  does  not  exist  in  a  pronounced  form  in  rodents, 
which  are  the  animal  models  most  frequently  used  for  studies  of 
the  molecular  mechanisms  underlying  synaptogenesis  [6], 

Common  marmosets  become  sexually  mature  earlier  (1.5  years 
of  age)  than  macaques  (sexual  maturation  at  approximately 
4.5  years  of  age  [8]).  It  is  possible  to  make  genetic  modifications 
in  the  marmoset  to  create  transgenic  animals  [9];  thus,  marmosets 
appear  to  be  a  convenient  model  for  examining  molecular  mecha¬ 
nisms  of  synapse  formation  in  the  primate  brain. 

In  a  previous  study,  we  examined  the  postnatal  development  of 
basal  dendritic  trees  and  spines  of  layer  III  pyramidal  cells  in  the 
primary  visual  sensory  cortex  (VI),  visual  association  cortex  (infe¬ 
rior  temporal  area,  IT),  and  the  vetrotalteral  prefrontal  cortex  (PFC) 
in  marmosets  [7]  and  reported  that  all  3  areas  show  overs  hoot-type 
synaptic  formation. 

In  this  study,  we  examined  gene  expression  of  3  cortical  areas  at 

5  postnatal  ages,  using  the  same  cortical  areas  and  timing  as  de¬ 
scribed  previously  [7],  on  a  custom-made  gene  chip  with  approxi¬ 
mately  23,000  marmoset  probes.  We  selected  statistically 
significant  differences  in  the  expression  of  genes  between  3  and 

6  months  (3  M  and  6  M)  of  age,  mostly  in  PFC,  IT,  and  VI,  for  iden¬ 
tifying  the  molecular  mechanisms  underlying  the  pruning  phase. 
Here,  we  report  and  discuss  the  selected  genes  that  participate  in 
axon  guidance  signaling,  which  is  the  most  plausible  (lowest  p  va¬ 
lue)  canonical  signaling  in  our  analysis  (Table  1 ).  We  also  discuss 
genes  of  the  microglia  and  immune  systems,  which  are  likely  to 


Canonical  pathways  most  likely  related  to  controlling  dendritic  spine  number,  based 
on  differentially  expressed  genes  between  3  and  6  months  of  age  in  the  3  brain 
regions. 


Red  letters  indicate  that  expression  at  3  months  was  stronger  than  at  6  months. 
Black  text  indicates  the  opposite. 

’  (Macrophage  system):  Genes  not  expressed  in  microglia  were  extracted  from  the 
original  canonical  pathway  list. 

f  Microglia-specific  genes  were  expressed  in  microglia  but  not  in  other 
macrophages. 


be  involved  in  synaptic  pruning  in  the  cerebral  cortex  and  in  plau¬ 
sible  canonical  signaling  (Table  1). 

2.  Materials  and  methods 

2.1.  Animals  and  tissue  preparation 

We  used  common  marmosets  (Callithrix  jacchus):  4  animals  on 
postnatal  Day  0  (0  D)  and  the  remaining  animals,  each  age  group 
consisting  of  3  animals  [2  M,  3  M,  6  M,  and  >1.5  years  old  (adult); 
Table  2],  All  experimental  and  animal  care  procedures  were  per¬ 
formed  in  accordance  with  the  N1H  guidelines  [11]  and  the  Guide 
for  Care  and  Use  of  Laboratory  Primates  published  by  the  National 
Institute  of  Neuroscience,  National  Center  of  Neurology  and  Psy¬ 
chiatry  (2005,  2008). 

2.2.  Tissue  processing 

The  animals  were  sedated  with  ketamine  hydrochloride  (Keta- 
lar,  25  mg/kg  i.m.;  Daiichi  Sankyo,  Tokyo,  Japan)  and  euthanized 
via  overdose  with  sodium  pentobarbital  (Somnopentyl,  100  mg/ 
kg  i.p.;  BCM  International,  Hillsborough,  NJ).  The  animals  were  per¬ 
fused  intracardially  with  0.1  M  potassium  phosphate-buffered  sal¬ 
ine  (pH  7.2).  Blocks  of  tissue  from  VI,  IT,  and  PFC  were  excised  as 
described  previously  [7],  and  immersed  in  RNAlater  (Life  Technol¬ 
ogies  Japan  Ltd.,  Tokyo,  Japan).  Total  RNA  was  extracted  using  the 
RNeasy  Mini  Kit  (Qiagen)  according  to  the  manufacturer’s  instruc¬ 
tions.  RNA  concentrations  in  the  samples  were  measured  using  the 
Amersham-Pharmacia  spectrophotometer  (Ultraspec  3100pro, 
Amersham  Pharmacia  Biotech,  Uppsala,  Sweden)  with  quantifica¬ 
tion  and  assessment  of  rRNA  integrity  using  a  2100  Bioanalyzer 
(Agilent,  Palo  Alto,  CA,  USA). 

2.3.  Generation  of  common  marmoset  cDNA  sequences  and 
construction  of  the  microarray 

The  DNA  sequences  for  probes  used  in  Marmo2  were  con¬ 
structed  as  follows.  The  human  RefSeq  sequences  were  BLAST 
searched  against  the  marmoset  6x  whole  genome  shotgun  preli¬ 
minary  assembly  sequence  (Ver.  2.0.2)  obtained  from  the  public 
domain  at  the  Genome  Center  of  Washington  University  School 
of  Medicine,  and  the  human  sequences  were  replaced  with  the 
marmoset  ones  by  means  of  the  BLAST  output.  After  excluding  re¬ 
peated  and  multiple  sequences,  19,117  human  orthologous  mar¬ 
moset  sequences  were  obtained.  Those  sequences  were 
combined  with  cDNA  sequences  used  for  the  design  of  the  Marmol 
microarray  [12]  and  newly  cloned  cDNAs  from  marmoset  ES  cells 
to  make  a  total  of  22,971  marmoset  sequences,  and  were  used 
for  designing  Affymetrix  GeneChip  probe  sets.  The  microarray 
was  constructed  according  to  the  GeneChip®  CustomExpress®  Ar¬ 
ray  Program  (Affymetrix,  Santa  Clara,  CA,  USA). 

2.4.  High-density  oligonucleotide  microarray  analysis 

Total  RNA  samples  were  prepared  using  the  RNeasy  Mini  Kit 
(Qiagen)  according  to  the  manufacturer’s  instructions,  and  concen¬ 
trations  were  measured  using  the  Nano  Drop  2000  spectrophotom¬ 
eter  (Thermo  Fisher  Scientific,  Wilmington,  DE).  rRNA  integrity  was 
assessed  using  a  2100  Bioanalyzer  (Agilent).  For  microarray 
probing,  reverse  transcription,  second-strand  synthesis,  and  probe 
generation  were  all  accomplished  with  the  3'IVT  Express  Kit 
(Affymetrix)  following  the  manufacturer’s  protocol.  In  brief,  from 
100  ng  of  total  RNA,  the  double-strand  cDNA  was  synthesized. 
From  the  double-strand  cDNA,  biotin-labeled  cRNA  was  prepared 
and  15  pg  of  labeled  cRNA  was  fragmented.  The  Marmo2a520631F 
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MA03 
MA01 
MAI  6 


Male 

Female 

Male 

Female 

Male 

Female 

Male 

Female 

Female 

Female 

Male 

Male 

Male 

Male 

Male 


e:  0  D,  0  days  of  age;  2  M,  2  months  of  age;  3  M,  3  months  of  age;  6  M,  6  months  of  age. 


Body  weight  (g) 
28.9 
32.0 


120.7 

96.2 


103.0 

142.7 

155.7 


298.5 

332.6 
370.2 


arrays  were  hybridized  as  described  in  the  Gene  Chip  Expression 
Analysis  Technical  Manual  (Affymetrix).  The  arrays  were  stained 
and  washed  using  the  GeneChip  Hybridization,  Wash,  and  Stain 
Kit  (Affymetrix),  while  fluorescence  was  intensified  using  the  anti¬ 
body  amplification  method  following  the  manufacturer’s  protocol. 
The  arrays  were  scanned  with  a  GeneChip  Scanner  3000  (Affyme¬ 
trix)  and  the  obtained  images  were  analyzed  using  Affymetrix 
microarray  suite  version  5  software  (MAS  5).  We  then  obtained 
tab-delimited  files  containing  data  on  the  relative  levels  of  expres¬ 
sion  of  transcripts  (Signal)  and  the  reliability  of  detection  (Detec¬ 
tion  Call).  The  derived  signal  values  were  globally  normalized 
and  the  control  for  all  probes  was  set  to  100. 


in  Supplemental  Table  1.  To  understand  the  biological  implications 
of  certain  expression  patterns,  these  data  were  analyzed  using  soft¬ 
ware  IPA  (Qiagen),  which  uses  curated  information  gathered  from 
scientific  literature.  We  further  confirmed  the  results  of  IPA  by  ver¬ 
ifying  original  research  articles. 

3.  Results  and  discussion 

Based  on  the  above  criteria,  571  genes  were  expressed  more 
strongly  at  3  M  than  at  6  M  and  1006  genes  were  expressed  less 
strongly  at  3  M  than  at  6  M  (Supplemental  Table  1 )  in  all  3  brain 
regions  examined. 


2.5.  Microarray  data  analysis 

To  examine  genes  specifically  expressed  in  different  cortical  re¬ 
gions,  average  signal  values  of  three  replicates  for  each  region  were 
calculated.  The  gene  expression  data  were  analyzed  using  Gene- 
Spring  GX  11  software  (Agilent).  The  signal  intensity  data  were 
preprocessed  using  the  RMA  algorithm  [13]  and  normalized  using 
Affymetrix  microarray  suite  version  5  software  (MAS5).  To  im¬ 
prove  the  false  discovery  rate  when  gene  probe  sets  were  turned 
on  or  off,  we  restricted  probe  sets  (>50%  present),  as  described 
by  McClintick  and  Edenberg  [14],  Differentially  expressed  genes 
among  the  3  areas  and  5  time  points  were  analyzed  using  2-way 
analysis  of  variance  (2-way  ANOVA;  p  value  computation:  asymp¬ 
totic;  multiple  testing  correction:  Benjamini  and  Hochberg  false 
discovery  rate;  cutoff  p  value:  1%).  In  2-way  ANOVA,  a  main  effect 
for  gene  expression  of  each  gene  was  either  only  age  or  only  brain 
region.  The  data  showing  a  main  effect  of  only  age  were  analyzed 
using  1-way  ANOVA  with  a  post  hoc  multiple  comparison  test  (Tu- 
key’s  HSD  method)  to  find  genes  differentially  expressed  between 
3  M  and  6  M  in  all  brain  regions  examined.  To  avoid  false  positive 
results  if  gene  distribution  was  potentially  asymmetrical,  we  used 
the  following  approach.  If  differentially  expressed  genes  were  fur¬ 
ther  filtered  so  that  the  absolute  value  of  FC  of  gene  expression  be¬ 
tween  3  M  and  6  M  is  greater  than  1.2,  this  method  efficiently 
damps  true  gene-regulatory  events  in  order  to  minimize  their  mis¬ 
leading  effect  on  the  normalization  process.  FC  was  calculated  as 
follows:  when  gene  expression  at  3  M  was  greater  than  that  at  6 
M,  the  former  was  divided  by  the  latter.  When  gene  expression 
at  6  M  was  greater  than  that  at  3  M,  FC  was  calculated  by  dividing 
the  former  by  the  latter  and  was  expressed  as  a  negative  value  (-). 
Thus,  negative  FC  values  indicate  that  gene  expression  increased 
from  3  M  to  6  M.  Information  on  the  selected  genes  is  presented 


3.1.  Axon  guidance  signaling 

In  the  canonical  pathway  analysis  of  IPA  using  the  list  above, 
axonal  guidance  signaling  (Supplemental  Table  2,  Table  1)  showed 
the  lowest  p  value.  Thus,  axonal  guidance  signaling  was  the  most 
plausible  pathway  based  on  genes  differentially  expressed 
between  3  M  and  6  M.  Although  the  list  of  canonical  pathways 
(Supplemental  Table  2)  according  to  our  data  included  axon  guid¬ 
ance-related  signaling  (Eph,  CXCR4,  IGF-1,  and  Wnt/p-catenin  sig¬ 
naling  pathways),  these  pathways  are  components  of  the  same 
system  of  axon  guidance  signaling.  Thus,  in  this  study,  the  genes 
that  we  discuss  also  include  these  subsystems  as  participants  in 
axon  guidance  signaling.  Table  3  shows  that  expression  of  genes 
belonging  to  axon  guidance  signaling  changes  from  3  M  to  6  M 
so  that  synapse/dendritic  spine  number  would  decrease  at  6  M  (Ta¬ 
ble  3,  [15-28]),  although  FC  is,  on  an  average,  rather  small.  These 
results  are  suggestive  of  some  role  of  axon  guidance  signaling  in 
pruning  approximately  at  6  M. 

3.2.  Immune  system 

Microphages  (part  of  microglia  in  the  brain)  and  proteins  that 
tag  substrates  in  synapses  to  be  eliminated  by  the  immune  system 
(e.g.,  complement  components)  play  an  important  role  in  synaptic 
pruning  [10],  Our  canonical  pathway  analysis,  as  mentioned  above, 
showed  that  the  macrophage  system  and  target-tagging  compo¬ 
nents  of  the  immune  system  were  the  pathways  that  were  most 
likely  to  undergo  changes  between  3  M  and  6  M  (macrophage  sys¬ 
tem,  hepatic  fibrosis/hepatic  stellate  cell  activation  as  the  3rd  most 
plausible  pathway;  dendritic  cell  maturation  pathway  as  the  84th 
plausible  pathway;  target-tagging  components  of  the  immune  sys¬ 
tem/complement  system  as  the  8th  most  plausible  pathway).  In 
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How  different  synapse  number-related  genes  belonging  to  axon  guidance  signaling  were  expressed  between  3  and  6  months,  and  expected  effects  of  different  expression  of  these 
genes  on  the  number  of  synapses/dendritic  spines. 


Pathway  within  axon  guidance  signaling  Possible  synapse  number-related  Fold  change  Expected  effect  of  gene  expression  change  on  Related 

genes  in  a  subpathway  (6  M/3  M)  synapse/dendritic  spine  number  at  6  M  compared  to  3  M  references 


Ephrin  signaling 


ROCK2 

RAC3 

EPHB2 

LIMK2 

RASA1 


1.27  Decrease 

1.21  Decrease 

1.47  Decrease 

1.20  Decrease 

-1.24  Decrease 


15 


19] 


SHANK2  signaling  RAPGEF7 

IGF-1  signaling  CDC42 

Tubulin-KLC-CRMP2  signaling  TUBA1C 

TUBB3 

TUBB2A/TUBB2B 

KLC1 

CXCR4  signaling  GNAZ 

Integrin  signaling  RASSF5 

Netrin/DCC  signaling  MKNK1 

Netrin/UNC/DCC  signaling  UNC5C 

NGF  receptor  signaling  PRKCG 

Wnt5a/FZD  signaling  PRKCG 

NFATC3 


1.21  Decrease 

1.34  Decrease 

1.21  Decrease 

1.21  Decrease 

1.43  Decrease 

1.23  Decrease 

1.20  Decrease 

-1.32  Decrease 

1.22  Decrease 

-1.25  Decrease 

1.30  Decrease 

1.30  Decrease 

1.23  Decrease 


this  section,  we  discuss  the  macrophage  system  and  the  target-tag¬ 
ging  components  of  the  immune  system  based  on  results  of  the 
canonical  pathway  analysis,  which  were  confirmed  in  the  brain 
only  via  the  synthesis  of  the  data  published  in  scientific  literature. 

3.2.1.  Microglia  (macrophage)  system 

3.2.1. 1.  AIF1.  Expression  of  AIF1  [29],  which  is  a  microglia-specific 
gene  that  reflects  the  activity  of  this  type  of  cells  [30],  decreased 
at  6  M  (Table  f). 

3.2.1. 2.  Microglia-modulating  genes.  Activity,  motility,  and  release 
of  cytokines  by  microglia  are  controlled  by  many  molecules  [10], 
Within  the  list  of  the  genes  differentially  expressed  between  3  M 
and  6  M  (Table  1),  CX3CR1,  CX3CR1,  TNFRSF19,  PDGFRA,  LEP,  LEPR, 
IGFBP3,  IGFBP4  TLR4,  EDNRB,  TIMP2,  COL1A2,  FN1,  1RF8,  PDGFRA, 
and  TREM2  were  downregulated  at  6  M  (FC:  1.22,  1.33,  1.26, 
1.31,  1.32,  1.36,  2.00,  1.34,  1.26,  1.37,  1.49,  1.92,  1.32  1.37,  1.31, 
and  1.24,  respectively).  TGRBR2,  which  is  an  inhibitor  of  microglial 
activity,  was  upregulated  (FC:  -1.23)  at  6  M.  Downregulation  or 
upregulation  of  the  above  genes  is  known  to  decrease  activity, 
motility,  and  the  release  of  cytokines  by  microglia  [31-43],  proba¬ 
bly  what  happened  at  6  M. 

3.2.2.  Target-tagging  components  of  the  immune  system 

3.2.2. 1.  Complement  system.  A  complement  protein,  C3b  partici¬ 
pates  in  opsonization  (enhances  phagocytosis  of  antigens  by 
microglia  in  the  brain)  [44],  C3-convertase  promotes  cleavage  of 
C3  into  C3a  and  C3b.  C1Q_B,  C1QC,  and  CFB,  which  are  involved  in 
the  synthesis  of  C3  convertase  [44],  were  upregulated  relatively 
strongly  (FC:  -1.52,  -2.24,  and  -2.05  respectively)  at  6  M. 

3.22.2.  Other  target-tagging  components  of  the  immune  system.  Both 
HLA-A  and  FCGR1B,  signals  for  the  elimination  of  synapses  by 
microglia  [45-47],  were  upregulated  (FC:  -1.28  and  -1.57,  respec¬ 
tively)  at  6  M. 

There  appeared  to  be  2  tendencies  of  gene  expression  related  to 
the  immune  system  between  3  M  and  6  M.  (1)  In  general,  microg¬ 
lial  activity  decreased  at  6  M.  This  result  suggests  that  although 
microglial  activity  is  thought  to  be  related  to  phagocytosis  of 
synapses  and  dendritic  spines,  microglial  activity  alone  could  not 
explain  how  pruning  was  accelerated  in  the  late  pruning  phase. 


(2)  Genes  (e.g.,  complement  components)  that  are  related  to  tag¬ 
ging  of  targets  of  microglial  synapse  pruning  tended  to  be  upregu¬ 
lated  relatively  strongly  at  6  M.  These  data  suggest  that  those 
genes  performed  important  functions  in  implementation  of  the 
pruning  phase  in  marmosets. 

4.  Conclusions 

In  summary,  the  axon  guidance  signaling  pathway  and  immune 
system  (a  combination  of  microglial  and  complement  components) 
probably  contributes  to  the  decrease  in  the  number  of  synapses/ 
spines  by  6  M  and  may  explain  the  overshoot-type  synaptic  forma¬ 
tion  in  primates.  Further  experiments  are  required  to  elucidate 
how  much  these  genes  and  subsystems  contribute  to  the  decrease 
in  the  spine  number  between  postnatal  3  M  and  6  M. 
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Appendix  A.  Supplementary  data 

Supplementary  data  associated  with  this  article  can  be  found,  in 
the  online  version,  at  http://dx.doi.Org/10.1016/j.bbrc.2014.01.024. 
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